Momentum distribution of ID Bose gases at the quasi-condensation crossover: 
theoretical and experimental investigation 
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We investigate the momentum distribution of weakly interacting ID Bose gases at thermal equi- 
librium both experimentally and theoretically. Momentum distribution of single ID Bose gases is 
measured using a focusing technique, whose resolution we improve via a levitation scheme. The mo- 
mentum distribution compares very well with quantum Monte Carlo calculations for the Lieb-Liniger 
model at finite temperature, allowing for an accurate thermometry of the gas that agrees with (and 
improves upon) the thermometry based on in situ density fluctuation measurements. The quasi- 
condensation crossover is investigated via two different experimental situations, each lying on either 
side of the crossover. Classical field theory is expected to correctly describe the quasi-condensation 
crossover of weakly interacting gases. We derive the condition of validity of the classical field theory, 
and find that, in typical experiments, interactions are too strong for classical field to be accurate. 
This is confirmed by a comparison between the classical field predictions and the exact quantum 
Monte Carlo results. 

PACS numbers: 03.75.Hh, 67.10.Ba 



I. INTRODUCTION 



Correlation functions are essential to describe many- 
body systems. In particular, the first-order corre- 
lation function, or equivalently its Fourier transform, 
the momentum distribution, is an important observ- 
able, since it witnesses various phenomena such as the 
Bose-Einstein condensation transition, the Berezinskii- 
Kosterlitz-Thouless transition [l[ or the Mott transi- 
tion 0]. In the case of one-dimensional (ID) homoge- 
neous gases, although one does not expect phase tran- 
sitions, the correlation functions contain valuable infor- 
mation about the system. While the equation of state of 
ID homogeneous Bose gases has been extensively stud- 
ied in various regimes |3H6| and compared with exact 
theories Q, the momentum distribution has remained 
largely unexplored. It was measured for quas i- ID quasi- 
condensates using Bragg spectroscopy |8l-fl0|. and more 
recently it was investigated using the focusing technique 
for quasi-lD Bose gases in the crossover from the ideal 
Bose gas to quasi-condensate regimes d Density 
ripples that appear at the near field of a freely ex- 
panding quasi- ID quasi-condensate have also been stud- 
ied fl2l . [l3l | . Although related to the second-order cor- 
relation function, they also provide information about 
the first-order correlation function within the quasi- 
condensate theory. The momentum distribution of an 
array of strongly correlated ID Bose gases [l4| consti- 
tutes the only measurement on a truly ID system to 
our knowledge. From a theoretical point of view, the 
momentum distribution of ID Bose gases with repulsive 
contact interactions is not known exactly over the entire 
phase diagram, but a few results have been established: 
the mean kinetic energy can be extracted from the ex- 
act Yang- Yang thermodynamics [ll| , and the short-range 



correlations are responsible for a momentum tail scaling 
as i/p 4 [IM3. The momentum distribution is known in 
the asymptotic regimes such as the ideal Bose gas and the 
quasi-condensate regimes, and a classical field approxi- 
mation has been prop osed to account for the crossover 
between them [U GlHll. 

In this work we investigate the momentum distribution 
of weakly interacting, purely or almost purely ID Bose 
gases in various regimes around the quasi-condensation 
crossover. On the theory side, we find a criterion for 
the classical field theory without cut-off [HI, [l^ to be 
accurate at the quasi-condensation crossover. Interac- 
tions are however too strong in our experiment to fulfill 
this criterion, and we implement exact Quantum Monte 
Carlo (QMC) simulations for the equilibrium behavior of 
the Lieb-Liniger model. We find that the measured mo- 
mentum distribution agrees very well with QMC calcu- 
lations. Moreover, the temperature extracted from a fit 
of the momentum distribution to the QMC calculations 
is in agreement with that obtained from in situ density 
fluctuation measurements 0, [|| . 

This paper is organised as follows. In Sec. [II] we re- 
call the classical field approximation and investigate the 
conditions under which it is valid to describe the quasi- 
condensation crossover. We show that, in most experi- 
mental situations, the gas is too strongly interacting for 
this description to be accurate. We then present our 
QMC calculations, and compare them with the classical 
field prediction in the asymptotic limit. In Sec. lIIIl we de- 
scribe the experimental setup and the focusing technique 
we use to measure momentum distribution. We discuss 
how we improve the resolution of this technique by means 
of a levitation scheme. In Sec. IIV1 we present our ex- 
perimental results that span both the quasi-condensate 
regime (purely ID) and the degenerate ideal Bose gas 
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regime (almost ID), and compare them to QMC results. 
We conclude and discuss the prospective work in Sec. [V] 



II. THEORETICAL PREDICTIONS 

A. Model Hamiltonian and exact solution 

We consider a ID homogeneous Bose gas with repul- 
sive contact interactions, described by the Lieb-Liniger 
Hamiltonian, whose grand-canonical expression is 
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dz 
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(1) 



where z is the position, Vl/ = ^(z) is the field operator, 
g is the coupling constant, m is the mass of one boson 
and \x is the chemical potential. The thermal equilib- 
rium state of the system is completely characterised by 
the density p and the temperature T, or by the dimen- 
sionless interaction parameter, 7 = and the reduced 

temperature, t 



2h i k B T 



The phase diagram Q in the 
parameter space (7, t) is given in Fig. [TJ where lines 
should be understood as crossovers. In this article, we 
focus on the quasi-condensation crossover, which occurs, 
for 7 < 1, around the line 7 co = t~ 2 / 3 (solid line). The 
Lieb-Liniger model has the remarkable property of be- 
ing integrable, giving access to exact results valid over 
the entire phase diagram. For example, the equation of 
state is known exactly through the Yang- Yang theory 0] . 
Yet, the calculation of the first order correlation function 
gW(z) = (* t (z)*(0)) from the Bethe-Ansatz solution 
for finite 7 and at finite temperature is still a subject of 
active research. Approximate analytic theories for the 
correlation functions exist in the limiting cases of ideal 
Bose gas regime (where the ideal Bose gas theory applies) 
and quasi-condensate regime (where the Bogoliubov the- 
ory applies [17jj). However, the behavior of g^\z) at the 
crossover between those two regimes lacks a quantitative 
description. 



B. Classical field theory 

The full quantum many-body problem of Eq. ([I]) is 
notoriously complex. However, since thermal fluctua- 
tions are expected to dominate at the quasi-condensation 
crossover, it seems appropriate to simplify the problem 
using a classical field approximation, where the quantum 
field operators ^5>(z) and W(z) are replaced by the com- 
plex fields and ^*(z) [H, GJ|. 

This strategy has 

been pursued in (ill . I20L [22] . Such an approach is partic- 
ularly attractive for ID homogeneous gases: contrary to 
higher dimensions, no ultraviolet divergence occurs when 
calculating g lyl \z), so that no energy cut-off is required 
for this particular quantity. This problem can then be 
mapped onto a time-independent 2D Schrodinger prob- 
lem, describing a single particle evolving in imaginary 
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FIG. 1: Phase diagram of the ID Bose gas in the param- 
eter space (7, i). The ideal Bose gas to quasi-condensation 
crossover occurs for 7 ~ 7 co = t~ 2 ^ z (solid line). Within 
the ideal Bose gas regime (light gray area), the crossover 
from the classical to the degenerate sub-regime occurs for 
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"id = t 



-1/2 



(dashed line). The region 7,( > 1 (dark 



gray area) is the fermionised regime. The segments corre- 
spond to the regions explored by the two sets of experimental 
data (see Sec. IIV[I . Note that the effect of quantum fluctua- 
tions on the function in the quasi-condensate regime are 
noticeable only for extremely low temperature [21| . which is 
experimentally irrelevant. 



time [18| , thus leading to fast and easy-to- implement cal- 
culations. 



The classical field model is parameterised by a sin- 



gle dimensionless variable v 
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vided that lengths are rescaled by the correlation length 
l c = h 2 pj (m/c^T) [l8l. [23(1 . Simple analytic formulas are 
found in the two limits of the ideal Bose gas regime 
(y 1) and the quasi-condensate regime (y <C 1): 
the momentum distribution is Lorentzian in both lim- 
its, with its full width at half maximum (FWHM) being 
Ap = 2h/l c in the ideal gas limit and Ap = h/l c in the 
quasi-condensate limit. One recovers here the expres- 
sions derived from the highly degenerate ideal Bose gas 
model and the high-temperature Bogoliubov theory |17| 
respectively. Results of numerical calculations spanning 
the crossover between those two regimes are reported in 
Fig. [21 where one assumes a fixed value t = 1000 and 
parameterises the system with 7 instead of v. We define 
the crossover region as the domain for which Ap (shown 
in Fig.rSfa) as the dashed line), differs by more than 10% 
from both asymptotic limits (h/l c and 2h/l c ), and we 
find that the crossover extends over about one order of 
magnitude in u, or equivalently in 7 at fixed t. Looking 
at the shape of the momentum distribution, we recover 
the Lorentzian distribution in both asymptotic regimes 
and, for any value of v, we find momentum tails decreas- 
ing as 1 /p 2 , as expected from the analog single-particle 
problem [24| . 
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C. Validity of the classical field theory at the 
quasi-condensation crossover 

The classical field theory is valid only if the population 
of the relevant modes is large. Within the ideal Bose 
gas regime, this requires p p& = sfrnksT /h, i.e. a 
highly degenerate gas. This condition writes equivalently 
7 *C 7d where 7^ = t~ x / 2 is the interaction parameter at 
degeneracy. Thus, the classical field theory only correctly 
describes the quasi-condensation crossover (which occurs 
at 7 ~ 7 co ) provided that 7 co <C 7<j. The last condition 
translates into i 1 / 6 3> 1. Since crossovers span typically 
about one order of magnitude in 7, one requires that 
t 10 6 - This value is very difficult to achieve experi- 
mentally on cold-atoms experiments while maintaining a 
temperature sufficiently low to ensure the ID condition, 
ksT <C hiuj_, unless extremely weak atomic interactions 
are reached, using for instance a Feshbach resonance (25[. 

In the experiment presented here, we have t < 1000. In 
this case, according to the above argument, the classical 
field approach is expected to be inaccurate at the quasi- 
condensation crossover. This can be seen in Fig. [Ha), 
which shows the FWHM of the momentum distribution 
as a function of 7 according to the ideal Bose gas theory 
and the classical field prediction for t = 1000. The ideal 
Bose gas theory is parameterised by \ — "f/jd — Pd/p, 
which quantifies the level of degeneracy: the FWHM of 
the momentum distribution within this theory goes from 
the Maxwell-Boltzmann prediction for x ^ 1 to 2h/l c for 
X« 1, and this crossover spans more than one order of 
magnitude in x, or equivalently in 7 at fixed t. The classi- 
cal field theory correctly describes the quasi-condensation 
crossover only if it shares a common plateau with the 
ideal Bose gas theory at Ap = 2H/l c in the degenerate 
ideal Bose gas regime, where -f co <C 7 <C 7<j. As seen 
in Fig. EJa), this is not the case at all for t = 1000: the 
highly degenerate ideal Bose gas regime is not very well 
identified for such a small value of t as far as the momen- 
tum distribution is concerned, so that the classical field 
theory does not correctly describe the quasi-condensation 
crossover. 



D. Quantum Monte Carlo approach 

For t < 10 6 , one needs therefore a more accurate treat- 
ment than the above-mentioned classical field approxi- 
mation. We make use of QMC simulations to provide 
numerically exact results for the Lieb-Liniger model at 
finite temperature. Bosons in continuous space can be 
efficiently simulated making use of Monte Carlo tech- 
niques. Variational and diffusion Monte Carlo results for 
the first-order correlation function at T = are reported 
in Ref. [26[. At finite temperature, one can in princi- 
ple make use of path-integral Monte Carlo simulations 
[27l . , although we are not aware of explicit calcula- 
tions for the Lieb-Liniger model. This latter technique 
requires the discretization of the imaginary-time dimen- 



sion (Trotter approximation), and its influence on the 
simulation results has to be carefully removed. Here, we 
adopt a rather complementary approach, discretizing the 
spatial dimension instead. Indeed, bosonic lattice models 
lend themselves to efficient QMC simulations free of any 
Trotter approximation. By discretizing the field operator 



$W(z) « b { p/y/a, where a is the lattice spacing, bj, b 



are bosonic operators at site j, and z 
the Bose-Hubbard Hamiltonian 



ja, we obtain 



j 

-(p,-Vj)hj , (2) 

with hj = bjbj, and parameters J = h 2 /(2ma 2 ), U = 
g/a, p = p — 2J. We also include the presence of an 
external potential V(z), discretized to give Vj. The Lieb- 
Liniger parameters 7 and t read 
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and t = 



4k B T ( J 
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(3) 



where n = (hj) = pa is the lattice filling. 

We study the Bose-Hubbard approximation to the 
Lieb-Liniger Hamiltonian by ma king use of Stochastic Se- 
ries Expansion (SSE) QMC [2^, |30| . extensively used to 
investigate lattice bosons. The Lieb-Liniger model is cor- 
rectly recovered when the lattice spacing is much smaller 
than the correlation length in the system. In the classical 
ideal Bose gas regime, this condition amounts to require 
AdB/a = \/47rJ/ (fcsT) ^> 1. On the other hand, in the 
degenerate ideal Bose gas regime and quasi-condensate 
regime, we need l c /a 3> 1, implying 2J/(Ujt) ^> 1 and 
4J/(U^t) 3> 1 respectively. The first-order correlation 
and the momentum distribution are efficiently sampled 
with SSE-QMC during the directed-loop update [29U30j. 
The lattice simulation reproduces faithfully the momen- 
tum distribution of the Lieb-Liniger model for momen- 
tum p <C hit j a. 

Fig. [2] shows QMC calculations for a homogeneous 
gas at t = 1000, for values of 7 that span the quasi- 
condensation crossover. In Fig. [D(a) the FWHM of the 
momentum distribution calculated with QMC is com- 
pared to that of the classical field and of the ideal Bose 
gas models. The QMC results follow the ideal Bose gas 
prediction almost until the FWHM of the ideal Bose gas 
prediction crosses that of the classical field prediction, 
which occurs for 7 = j c , and then converges towards the 
classical field prediction. The disagreement to the re- 
spective theories never exceeds 20%. These results miti- 
gate the "failure" of the classical field approximation: as 
long as the FWHM is concerned, an approximate model 
where one uses the ideal Bose gas prediction for 7 > j c 
and the classical field prediction for 7 < j c , would give 
the correct predictions within 20% error. Investigation 
of the full momentum distribution is shown in Fig. fflb) 
for three different values of 7. On the ideal Bose gas side 
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FIG. 2: Momentum distribution of a homogeneous gas at a t = 1000. (a) FWHM according to the classical field theory (dashed 
line), ideal Bose gas theory (solid line), Maxwell-Boltzmann prediction (dash-dotted line) and QMC calculations (stars). The 
arrows show the interaction parameter at the quasi-condensation crossover (7 co = t~ 2 ^ 3 ) and at the degeneracy crossover 
(jd = t -1//2 ). (b) momentum distribution in log scale for three different value of 7. Shown are QMC results (solid lines), 
classical field predictions (dashed lines) and ideal Bose gas prediction (dotted lines). 



of the quasi-condensation crossover, for 7 = 5^/ co , the 
QMC momentum distribution follows closely the ideal 
Bose gas behaviour. The classical field result is recovered 
for 7 w 7 C o/3, while the momentum width is still about 
10% above the Bogoliubov prediction. For 7 = -f co , we 
find that the tails of the momentum distribution agree 
well with the ideal Bose gas prediction, while the width 
of the distribution is narrowed by about 20%. Finally, 
note that, although the highly degenerate ideal Bose gas 
prediction Ap = 2h/l c has not been reached before the 
gas undergoes the quasi-condensation transition, the mo- 
mentum distribution at 7 ~ 7 co is much narrower than 
the Maxwell-Boltzmann prediction (dash-dotted line in 
Fig. [2f a)) , so that the effect of degeneracy within the 
ideal Bose gas regime is already substantial. 



densities, when almost the whole cloud lies within the 
quasi-condensate regime, the LDA condition gives lo z <C 
ksT \/2mgj pjf? @. Both conditions are fulfilled in the 
data shown below, so that LDA is expected to be valid. 
Using QMC, we compute the momentum distribution of 
a harmonically trapped gas for different peak chemical 
potentials (i.e. different peak linear densities) and tem- 
peratures. To shorten the computational time, calcula- 
tions are performed for trapping frequencies five times 
larger than the experimental values, while still satisfying 
the LDA condition, as verified numerically. Linear in- 
terpolation between the calculated distribution permits 
the calculation of momentum distribution for any chem- 
ical potential [1q and temperature T, and enables us to 
perform fits to the measured distribution. 



E. Local density approximation 

While the previous results hold for homogeneous infi- 
nite systems, in the experiment described below, the gas 
is trapped longitudinally in a harmonic confinement of 
frequency lu z /(2tt). As long as the correlation length l c is 
much smaller than the cloud size L, however, a local den- 
sity approximation (LDA) is valid: at a given position z, 
the local value of an observable is that of a homogeneous 
gas at the chemical potential ^l(z) = /xq — m{<jj z z) 2 /2, and 
the global value for a trapped system can be obtained by 
adding the contributions of each slice of the gas, i.e., 

Otrap A») = 7= = / — OhomO - P) ■ (4) 

^2mujj Jo \i 

Therefore, within this approximation and for a given 
peak density, the normalized momentum distribution is 
independent of lo z , while the total atom number scales 
as 1/lo. For a gas whose peak linear density lies in the 
crossover from the degenerate ideal Bose gas regime to 
the quasi-condensate regime, the LDA condition l c -C L 

writes co z < (mg 2 k 2 B T 2 /h 5 ) 1/3 At much higher 



III. EXPERIMENTAL TECHNIQUES 

One-dimensional atomic clouds are realised in our ex- 
periment using an atom-chip setup (3l| . More specifi- 
cally, we load 8 x 10 4 Rb 87 atoms in the modulated guide 
described in The longitudinal harmonic potential has 
an oscillation frequency of about 8 Hz and is obtained by 
DC currents running through wires I3 and I4 (see Fig. [3]). 
The transverse potential is provided by three wires car- 
rying a 200 kHz sinusoidal current. The transverse oscil- 
lation frequency can be varied from a few hundred Hz 
to tens of kHz by changing the current amplitude. An 
important feature of this design is the independence be- 
tween the longitudinal and transverse confinements. Af- 
ter 1200 ms of radio-frequency evaporation, we let the 
cloud reach thermal equilibrium for 400 ms, after which 
no remaining breathing can be observed. To characterise 
the atomic cloud, we perform in situ density fluctuation 
measurements, analysing the noise in an ensemble of ab- 
sorption images as described in [||. 

To measure the momentum distribution, one can moni- 
tor the free expansion of the gas after all confining poten- 
tials have been turned off. The rapid transverse expan- 
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FIG. 3: (a) Wire configuration: currents Ji, /2, I3 and I4 re- 
alise the longitudinal trapping potential of frequency f z . The 
transverse confinement of frequency f± is ensured by the three 
wires carrying a current modulated at 200 kHz. (b,c) Focus- 
ing sequences. The focusing potential, harmonic to order 5 in 
z, is applied during tkick before the longitudinal potential is 
turned off. After a free evolution during the focusing time t, 
an absorption image is taken. Without levitation, see (b), the 
transverse confinement is removed during the whole focusing 
time, though a 1 ms ramp down to f± — 700 Hz is carried 
out to decrease the transverse expansion. In the presence of 
levitation, see (c), f± is kept at 700 Hz for a time ti, during 
which the cloud does not fall. Numerical values in (b) and (c) 
correspond to data of Fig. [5] and [4] respectively. 



sion amounts to an effective instantaneous suppression 
of the interactions, so that the longitudinal expansion 
reflects faithfully the initial momentum distribution in 
the trap. However, reaching the far field regime where 
the longitudinal density profile is homothetic to the mo- 
mentum distribution requires unrealistically large field of 
view and long expansion time for our experimental pa- 
rameters. We thus use the so called focusing technique, 
already applied to ID [Til. I32I |33| and 2D systems [3~4| . 
Its principle is described below. A strong longitudinal 
harmonic potential is applied for a very short time dur- 
ing which atoms do not have time to move but acquire 
a longitudinal momentum shift Sp = —Az, proportional 
to their distance z to the trap center. Then, the atomic 
interactions are suppressed by switching off the trans- 
verse confinement, and the longitudinal confinement is 
removed so that the cloud starts a free evolution. After 
a focusing time t = m/A, the longitudinal density distri- 
bution f(z) is homothetic to the initial momentum dis- 
tribution n(p) = (m / V) / '{jpr I 'm) [35j|. More precisely, the 
time sequence we use is drawn in Fig. [3] The longitudinal 
focusing potential is realized by a four-wire configuration 
that cancels out the non harmonic terms up to the 6th 
order, thus minimizing aberrations. It has an oscillation 
frequency Wkick/27r = 40 Hz and it is applied for about 
^kick = 0.6 ms. After the focusing pulse, we ramp the 
transverse frequency down to 700 Hz in 1 ms before com- 
pletely turning off the transverse trapping potentials and 
letting the cloud evolve freely. The 1 ms transverse ramp 
is used to limit the transverse velocity spread, so that the 
transverse final expansion of the cloud is reduced and the 



signal to noise ratio (SNR) is increased. We verified that 
this ramp is however quick enough to leave the longitudi- 
nal velocity distribution unchanged. After free evolution 
during the focusing time t ~ l/ w kick^kick ~ 27 ms, we 
take an absorption image of the cloud and extract the 
longitudinal density distribution, from which we deduce 
the initial momentum distribution. The focusing time r 
is adjusted by minimising the width of the density dis- 
tribution of cold samples. Moreover, the width of very 
cold clouds saturates at a value that gives us the mo- 
mentum resolution Ap/h = 0.2 fj,m~ . This resolution is 
significantly smaller than the width of the data in Fig. [5] 

On the other hand, a better resolution may be neces- 
sary for a momentum distribution that is much narrower. 
This can be achieved by a longer focusing time r. How- 
ever, the depth of field of our optical system limits the 
maximum free fall under gravity to about 3 mm, corre- 
sponding to a free fall time of about 25 ms. Nevertheless, 
larger focusing times can be obtained if a transverse po- 
tential is maintained during the focusing time, strong 
enough to overcome gravity (which acts in the trans- 
verse plane in our experiment), but weak enough so that 
atomic interactions are negligible. More precisely the lev- 
itation scheme we use is the following. After having per- 
formed the longitudinal focusing pulse and switched off 
the longitudinal potential, we ramp the transverse fre- 
quency down in 1 ms, similar to the procedure described 
before. The final transverse frequency is /; = 700 Hz, 
which is the lowest value that allows the atoms to remain 
transversally trapped. We then hold the cloud for a levi- 
tation time ti , at the end of which we switch off the trans- 
verse confinement and let the cloud fall freely for t f . The 
levitation does not affect the longitudinal free evolution 
as long as the effect of interactions stays negligible. We 
have checked that this is indeed the case for the explored 
parameter range by measuring the longitudinal density 
profile at focus f(z) of identical clouds, either without 
levitation and a focusing time r = tf = 24.7 ms, or with 
the above levitation scheme and r = t; + tf = 43.3 ms, 
where t\ = 20.6 ms and tf = 22.7 ms. The functions 
/(mz/r), shown on Fig. [4] for a cloud of 7000 atoms at 
T w 95 nK, give indeed the same momentum distribu- 
tion. The validation of this technique paves the way to- 
wards high resolution measurements of momentum dis- 
tribution of ID gases. 



IV. EXPERIMENTAL RESULTS 



Fig. [5] shows experimental results referring to two dif- 
ferent parameter sets (particle number, temperature and 
trapping potential). In both cases, the momentum dis- 
tribution is obtained without the levitation scheme, with 
a focusing time of 27 ms. 
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FIG. 4: Momentum distribution of a gas of about 7000 atoms 
at thermal equilibrium, in a 7.5 Hz longitudinal trap, and 
4.6 kHz transverse trap. The solid distribution is obtained 
without levitation and with a free fall of 24.7 ms. The circles 
are the momentum distribution of the same sample, but with 
a levitation time of 20.6 ms followed by a free fall of 22.7 ms. 



A. On the quasi-condensate side of the crossover 

Data in Fig. HJA) correspond to a cloud ini- 
tially trapped with a transverse oscillation frequency 
wj_/(27r) = 6.4 kHz and a longitudinal frequency of 
8.3 Hz, measured using parametric heating and dipole os- 
cillations respectively. Since the cloud temperature and 
chemical potential are much smaller than the transverse 
energy gap, the system is well described by a ID gas with 
a coupling constant g = 2huj±a s , where a s = 5.3 nm 
is the 3D scattering length [11]. In Fig. (A, 1-2), we 
show the measured momentum distribution (points) with 
a QMC fit (solid lines) that yields T = 72 nK, corre- 
sponding to t = 76. Note that the shape is reproduced 
up to a few per cent. We report in Fig. [T]the segment 
[(7m j t), (7M) t)} corresponding to the domain explored by 
the data: 7 m corresponds to the peak density, and the 
segment length is such that 80% of the atoms lie in a 
region of density larger than that associated to 7m. We 
find that 7 m ~ 7co/2, where 7 co = i~ 2 / 3 is the interaction 
parameter at the quasi-condensation crossover. Thus the 
cloud lies quite deeply in the quasi-condensate regime. 
However, the measured momentum distribution is sub- 
stantially broader than that predicted using the quasi- 
condensate Lorentzian n(p) cx l/(p 2 + (h/l c ) 2 ) together 
with the LDA, shown as the dashed line in Fig. |3{A,1). 
On the other hand, for such a value of j m , the momen- 
tum distribution shows slowly decreasing tails, compati- 
ble with a Lorentzian behavior within a few percent in a 
wide range of momenta. This is seen in Fig.[5f A,2) where 
the dashed-dotted line shows a 1/p 2 decrease. For our 
SNR of about 80, data are compatible with a Lorentzian 
behavior at large p. 

As shown in [HQ, an independent and reliable ther- 
mometry can be carried out using in situ density fluc- 
tuations, provided the pixel size A is both much larger 
than the correlation length of density fluctuations, and 
much smaller than the typical length scale of the vari- 
ation of the mean density profile : the atom number 
fluctuations are then given by the thermodynamic re- 
lation (5N 2 ) = k B TAdp/dp, where p(p,T) is the equa- 
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FIG. 5: Data for atoms confined in a harmonic trap of trans- 
verse oscillation frequency of (A) 6.4 kHz and (B) 2.1 kHz, 
and a longitudinal oscillation frequency of (A) 8.3 Hz and 
(B) 7.6 Hz. The pixel size for in situ data is A = 2.7 pm. 
1),2) momentum distribution in linear and log scales. QMC 
fits are shown in solid lines, yielding (A) T = 72 nK and 
(B) T = 84 nK. The dashed line in (A,l) is the prediction 
using the LDA and the Lorentzian momentum distribution 
in the quasi-condensate limit. In 2), the dashed-dotted lines 
give the 1/p 2 behavior, while the momentum associated with 
the healing length £ = h/^/mpg indicates that sample (A) 
is more in the quasi-condensate, with Ap <g h/£, and sam- 
ple (B) is more in the ideal Bose gas, with Ap ^> h/£. In 
addition, the dotted lines indicate the contribution from the 
transverse excited states, which is insignificant for data (A) 
and hence neglected in our analysis. 3) in situ density fluc- 
tuations, with Yang- Yang prediction in solid lines, evaluated 
at the above-mentioned temperatures. The gray dashed and 
dotted lines shows Yang- Yang predictions with 30% tempera- 
ture deviations. The dash-dotted lines indicate the poissonian 
fluctuation level. 4) in situ profile, with the QMC profile in 
dashed lines. The dotted lines again show the contribution 
from transverse excited states. We also plot in solid lines the 
density profile fitted with Yang- Yang calculation, which yields 
(A) T = 111 nK and (B) T = 76 nK. 



tion of state of homogeneous ID gases, known exactly 
using the Yang- Yang calculation [6j, and the derivative 
is computed for [i such that p(fi,T) = (N) / A is the lo- 
cal density. The above criterion is well fulfilled in our 
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experiment and the measured density fluctuations are 
shown in Fig. [5jA,3). For linear densities larger than 
25 atoms per pixel, fluctuations are almost independent 
of the density, which confirms that the center part of the 
gas lies within the quasi-condensate regime. Indeed, in 
this regime [i ~ gp, so that the thermodynamic relation 
reduces to (SN 2 ) = AksT/g, which does not depend on 
(N). The fluctuations computed using the thermody- 
namic relation and the Yang- Yang equation of state, for 
the temperature T = 72 nK (obtained by fitting the mo- 
mentum distribution with the QMC results) are shown 
as solid line. The data are in agreement with this predic- 
tion. With our SNR, the density fluctuations based ther- 
mometry has a precision of about 20% (see gray curves in 
Fig.[5fA,3)), less precise than the thermometry obtained 
by fitting QMC calculations to the momentum distribu- 
tion. 

Finally, we compare in Fig. A, 4) the measured in situ 
density profile to the QMC density profile (within LDA) 
for T — 72 nK (dashed line). We find a good agreement 
for most of the profile, although the measured data show 
higher wings. Since the central part of the cloud lies 
deep in the the quasi-condensate regime, a large part of 
the cloud follows the Thomas-Fermi profile, and all the 
information on the temperature lies in the small wings. 
This renders Yang- Yang thermometry based on density 
profile less precise and extremely sensitive to the shape 
of the wings. Here, a Yang- Yang fit to the density pro- 
file gives T = 110 nK, a value 40% higher than that 
extracted from the QMC fit, incompatible with the mea- 
sured momentum distribution or the density fluctuation 
measurements. This discrepancy, and thus the presence 
of the inflated wings, may come from the anharmonicity 
of the potential due to its residual roughness (37|. Al- 
ternatively, it may also indicate a lack of perfect thermal 
equilibrium. 



B. On the ideal Bose gas side of the crossover 

While the previous results probe mainly the quasi- 
condensate regime, we also probe the ideal Bose gas 
side of the quasi-condensation crossover, i.e. data with 
7m > 7co = t~ 2 / 3 . The ideal Bose gas regime show a very 
different behaviour from the trivial Maxwell-Boltzmann 
prediction only for a large t parameter, for which the 
quasi-condensation crossover occurs for an already highly 
degenerate gas. If one wishes to preserve the ID condi- 
tion ksT <C hujj_, large t parameters can be accessed 
only by decreasing the transverse confinement (25) . We 
thus reduced the transverse confinement to 2.1 kHz. 
Data are shown in Fig. EJB). No saturation of the den- 
sity fluctuations is seen on the in situ fluctuation mea- 
surements, which indicates that the gas does not lie in 
the quasi-condensate regime. Fluctuations however rise 
well above the poissonian level (shown as dashed-dotted 
line on Fig. [5p3,3)) so that the gas is highly degener- 
ate. Contrary to the data of Fig. [H A), we now have 



fcsT '/ '(Hu>±) = 0.8 which is of the order of unity, so that 
the transverse excited states contribute to the measured 
fluctuations and momentum distribution. We take into 
account the population of transverse excited states, as- 
suming that they behave as independent ideal ID Bose 
gases, while the ground state is treated with QMC cal- 
culation. This modified QMC model has been used with 
success to describe density profiles [|| and density fluctu- 
ations Q , and has been applied to predict the rms width 
of momentum distribution 

|llj|. A fit using the above 
strategy reproduces well the measured momentum dis- 
tribution (see Fig. [S£B,1)), and yields the temperature 
T = 84 nK, corresponding to t = 840. The segment 
in the phase space (7, t) explored by the data is shown 
in Fig. [TJ the peak linear density is close to the quasi- 
condensation crossover density and most of the cloud lies 
in the degenerate ideal Bose gas regime. The contribu- 
tion of the excited states to the momentum distribution 
is seen in Fig. |5p3,l): it is only 10% in the center but it 
rises to almost 50% in the wings around \p\/H ~ 5 /xm _1 . 
Fig. [5jB,2) shows that the theoretical momentum dis- 
tribution decreases faster than 1/p 2 at large momenta. 
This is mainly due to the contribution of the excited 
states, which have approximately Gaussian momentum 
tails. However, with our SNR of about 50, no devia- 
tion from a 1/p 2 behavior can be identified. We report 
in Fig. [51(B, 2) the fluctuations expected for the tempera- 
ture T = 84 nK (obtained from the fit of the momentum 
distribution with QMC). They are in agreement with the 
measured fluctuations. Note however that, for these pa- 
rameters, the uncertainty of thermometry based on fluc- 
tuations is about 30% (see gray lines in Fig.[5jB,3)). The 
density profile expected for T = 84 nK is not far from the 
measured one (see Fig. |SjB,4)). In contrast to case (A), 
Yang- Yang thermometry based on the profile is less sen- 
sitive to tiny modifications of the wings, since the profile 
of the high-density regions is also affected by the tem- 
perature. Thus, a thermometry based on the profile is 
expected to be quite precise and the potential roughness 
is expected to have a smaller effect. A Yang- Yang fit to 
the experimental profile yields a temperature of 76 nK, 
different only by 10% from the temperature deduced from 
the momentum distribution. 



V. CONCLUSION AND PROSPECT 

We have shown that for t < 10 6 , a simple ID classi- 
cal field approximation without cut-off fails to describe 
the momentum distribution of weakly interacting ID 
Bose gases in the crossover from ideal Bose gas to quasi- 
condensate. Thus, we performed QMC calculations to in- 
vestigate the crossover. Experimentally, we measure mo- 
mentum distribution with the focusing technique, which 
is improve by a levitation scheme for enhanced resolu- 
tion. We show that the temperatures deduced from the 
momentum distribution is in agreement with an indepen- 
dent fluctuation-based thermometry. We find that the 
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naive Bogoliubov result is not appropriate for parame- 
ters that are close to the ones relevant to this article, even 
when the cloud lies quite deeply in the quasi-condensate 
regime. In it has been proposed to use the mean 
kinetic energy deduced from measured momentum dis- 
tribution in order to extract the temperature using the 
modified Yang- Yang calculation. For the data presented 
in this article, with a SNR of about 50, the measured mo- 
mentum distribution shows tails compatible with a 1/p 2 
behavior both in the quasi-condensate regime and in the 
degenerate ideal Bose gas regime. Therefore, extraction 
of the mean kinetic energy from the data is impossible. 
We argue that this technique was applicable in jTTj] be- 
cause of the more 3D nature of the gas: for quasi ID 
discussed in section IIVB[ the contribution of 
the transverse excited states leads to fast decaying tails, 
which enable the extraction of the mean kinetic energy 
per particle. We also remark that neither our measure- 
ments nor the QMC calculations we perform display the 
asymptotic behaviour of 1/p , which only appears at even 
larger momenta for our parameters. On the other hand, 
we believe that this behaviour could be more readily dis- 
cernable in the strongly interacting regime, where the 
1/p tail should contain a larger proportion of atoms. 

This work opens many perspectives for the study of 
ID gases. First of all, we show that momentum distri- 
bution measurements provide a precise thermometry for 
ID Bose gases. While in this article we use QMC cal- 
culations to fit the temperature of experimental data, a 
rough estimate of the cloud temperature could be per- 
formed with a lower numerical cost. Indeed, the compar- 
ison with the exact QMC momentum distribution show 



that a model which combines the ideal Bose gas theory 
and the classical field approximation gives the correct 
width within 20% precision. Second, momentum distri- 
bution measurements are essential to characterize more 
complex systems. For instance, in the presence of a lat- 
tice, it enables the investigation of the correlation prop- 
erties at the Mott and /or pinning transitions. Contrary 
to Bragg spectroscopy [lfj, Il4| , where a single momentum 
component is probed at each shot, the focusing method as 
well as time of flight method gives access to the whole dis- 
tribution at the same time, which allows for noise correla- 
tion measurements in momentum space (23l [38|. Finally, 
recording the time evolution of the momentum distribu- 
tion is essential to monitor out of equilibrium dynamics 
and address the questions about thermalisation in ID. 
For example, the dynamics resulting from a quench of 
the ID coupling constant can be investigated. The mea- 
sure of the momentum distribution of impurities would 
also permit the study of impurity dynamics, which is cur- 
rently under intense investigation [39j. 
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